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Abstract 

The formula for probability density functions (PDFs) has been extended to include PDF for 
energy dissipation rates in addition to other PDFs such as for velocity fluctuations, velocity deriva- 
tives, fluid particle accelerations, energy transfer rates, etc, and it is shown that the formula actually 
explains various PDFs extracted from direct numerical simulations and experiments performed in a 
wind tunnel. It is also shown that the formula with appropriate zooming increment corresponding 
to experimental situation gives a new route to obtain the scaling exponents of velocity structure 
function, including intermittency exponent, out of PDFs of velocity fluctuations. 

PACS numbers: 47.27.-i, 47.53. +n, 47.52. +j, 05.90.+m 

Keywords: intermittency, multifractal, singularities, scaling exponents, probability density function 



*T. Arimitsu: tarimtsu@sakura.cc.tsukuba.ac.jp 



1 



I. INTRODUCTION 



The quest for an essence of intermittency, i.e., a fundamental process in turbulence, has a 
long history but still an unsolved problem in physics for more than 120 years since about 1880 
when the systematic experiments of turbulence was started by Reynolds. The theoretical 



research on the subject in fully developed turbulence (we simply call it turbulence in 



lowing unless it is confusing) starts with Kolmogorov's dimensional analysis (K41) [l] based 
on the assumption of the self-similarity of fluctuating velocity field in the inertial range. 
After Landau's criticism against K41 in about 1944 and the preliminary research by Heisen- 
berg [2] , the quest develops, mainly, into two directions. One is the dynamical approach; the 
other is the ensemble approach [3j. Within the dynamical approach one treats the stochastic 
Navier-Stokes equation by perturbational methods whereas within the ensemble approach 
one performs statistical mechanical analysis of turbulence under the assumption that eddies 
make up energy cascade. It has been, gradually, revealed that, among 

the ensemble methods, a new theoretical framework named multifractal probability density 
function analysis (MPDFA) and A&A model within the framework can analyze in a high 
precision the data extracted out from the recent experiments and simulations conducted 
with higher accuracy. 

Several quantities such as velocity derivatives, fluid particle accelerations and energy 
dissipation rates have some singularities due to the invariance of the Navier-Stokes equation 
for velocity field u(x, t) in high Reynolds numbers under the scale transformation: 

x^x' = \x, u^u' = X a/3 u, t —>t' = \ 1 ~ a/3 t, p^p' = \ 2a / 3 p (1) 

with arbitrary real number a Here, p represents pressure. The MPDFA is a statisti- 
cal mechanical theory of an ensemble providing analytical formulae for various probability 
density functions (PDFs) applicable to intermittent systems. It was constructed on the as- 
sumption that the strengths of the singularities distribute themselves in a multifractal way in 
real physical space. This distribution of singularities determines the tail part of PDFs. The 
parameters appeared in the theory are determined, uniquely, by the intermittency exponent 
/j, that represents the strength of intermittency. 

On the other hand, observed PDFs should include the effect resulted from the term in 
the Navier-Stokes equation that violates the invariance under the scale transformation (the 
dissipative term). There has been, however, no ensemble theory of turbulence taking this 



the fol- 



effect into account, and the situation remained at the stage where almost all the theories 
are just trying to explain observed scaling exponents ( m of the mth order velocity structure 
function, i.e., the mth moment of velocity fluctuations. The MPDFA counts this effect as 
something determining the central part of PDFs narrower than its standard deviation. We 
are assuming that the fat-tail part, which the PDFs of intermittent systems took on, is 
determined by the global characteristics of the system, and that the central part of PDFs is 
a reflection of the local nature of constituting eddies. 



II. FORMULA FOR A&A MODEL WITHIN MPDFA 



The scaling exponents of the velocity structure function within A&A model is given by 

DO 



( m = a m/3 - 2Xm 2 / [9 (l + sjc^)] [l - log 2 (l + sjc^)] /(l - q) (2) 

with C m = 1 + 2m 2 (l — q)X\a2. The parameters «o, X and q are introduced through the 
Tsallis-type distribution function 

P {n \a)da oc {l - [{a - a ) / /ia] 2 }™ 7 ^ da (3) 

with {Aa) 2 = 2X /[(l — g)ln2] adopted within A&A model as the probability to find a 
singularity specified by a within the range a ~ a + da for large n, and are determined 
self-consistently as functions of \i through the relation fi = 2 — £g. This determines the 
distribution of the arbitrary real number a appeared in the scale transformation |T|). 

When the scaling exponents are given by numerical or ordinary experiments, we analyze 
them with the formula (J2J) to determine the value fi. When they are not given experimentally, 
we have another route to determine its value with the help of the observed PDFs. The latter 
new route is provided first in this paper within MPDFA. 

We list here a unified formula for PDFs H^\x n ) within A&A model of MPDFA. The 
contribution of singularities to PDF is taken into account by 

n^(\x n \)d(\x n \) P ( - n \a)da (4) 

with the transformation of variables 

\x n \ = 6x n /6x = 8f\ 6 n = £ n /£ = 2"", 6x n = \x(m + Q - x{m)\. (5) 



Here, x represents an observable quantity such as a component of fluid velocity field u, 
pressure p, etc., and n does a number of multifractal steps whose increment An gives the 
zooming increment that should correspond to the process how experimentalists extracted 
PDFs by changing the consecutive distances r = £ n between two observing points, say r' 
and r, i.e., with an appropriate /i the correct zooming increment An = n' — n is provided 
by An = — log 2 (r'/r) with n = —log 2 {r/f]) +log 2 {£o/f]) where n is the Kolmogolov length. 
Note that £q is a reference length that is not necessarily equal to the integral length in 
general 

The tail part (|£ n | > £*) of the PDF for variable £ n defined in the ranges (—00, 00) and 
(0, 00) is given by 



n ^tl(£nK" = tif,s{x n )dx n OC |^ 



(n) 



in 



1 — g (3 In 1^/^,0 



n 



n/{l-q) 



din- 



(6) 



20 2 X|m<5 n | 

The center part (|£ n | < £*) of the PDF for variable £ n defined in the range (—00, 00) is given 



by 



n&(e») oc <i - (i - 4 



+ 3/'(a*) 
20 



(7) 



whereas for variable £ n defined in the range (0, 00) by 



n 



in 



,cr(Cn) OC I — 



+ 3/V) 
20 



1 - 



1/(1-9') 



(8) 



The variable ^ n is the scaled variable related to an observed variable 5x n through the relation 
in = bx n j (Sxl) 1 ^ 2 . The tail part and the center part of PDFs are connected at £* under the 
conditions that they have a common value and a common log-slope. The point £* has the 
characteristics that the dependence of PDF on n is minimum for large n. We will see that 
i* ~ 1 through the analyses of experiments. The tail part (JOJ) is determined by (JBJ) with 
the translation of variable given by the first equation of (J5|i. Note that the formulae (jHJ) 
and (J7J) or (JSJ) are unified in the sense that it provides the PDFs of velocity fluctuations and 
of velocity derivatives with = 1, the PDFs of pressure fluctuations and of fluid particle 
accelerations with = 2, and the PDFs of energy transfer rates and of energy dissipation 
rates with = 3. Note also that the energy dissipation rate is a variable taking only positive 
real values, whereas the others are variables taking both negative and positive real values. 
The PDFs for the latter variables are given by (jOJ and (JJJ) with three parameters q, n and 
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q' . On the other hand, the PDF for the former variable is given by © and (JHJ) with four 
parameters q, n, q' and 9. 

These parameters are determined by the following procedure through a series of observed 
PDFs obtained by changing the distances of two measuring points, i.e., 1) Start with a trial 
value fi (and also with trial values q' and/or 9) to fit one of the observed PDFs with the 
tail part PDF given by ©. Note that the values of parameters «o, X and q are determined 
as functions of /i, self-consistently. 2) Once one has an appropriate fi value, other observed 
PDFs in the series can be fit with the correct increment An. 3) After getting the parameters 

(or equivalently q) and n with trial values q' and/or 9, one can adjust better values for 
q' and/or 9 by fitting the center part of the observed PDFs with the formula (J7J) or (JHJ). 4) 
Repeating the above process 1) ~ 3), one can obtain the best set of parameters. 

With the PDFs © and (J7J) or (jHJ), the mth moment of the quantity \x n \ is given by 

<kr)^7jl+(l-^) a^e" (9) 

where k = 2 for the variables with the range (— oo, oo) and k — 1 for (0, oo), a 3 ^ m = 
[2/y/C^(l + \JcZ)] l/2 and 

itl = J o °° dx n \x n r[uf\x n ) - <>(*„)]• (10) 

Since the scaling exponents Q is related to the generalized dimension Dq by |3] 

( 3 g = 1 - (1 - q)Dg, (-oo<g<oo) (11) 

one can determine the generalized dimension of the system through Q . 



III. ANALYSES OF SIMULATIONS AND EXPERIMENTS 

In Fig. 4 of and in Fig. 1 of 0, the PDFs of transverse velocity fluctuations and of 
transverse velocity derivatives measured in the DNS on 2048 3 mesh size ^| at = 380 are 
analyzed by the present theoretical PDFs (JHJ and (JZJ) with <fi = 1 for velocity fluctuations 
and for velocity derivatives, and are plotted on (a) log and (b) linear scales. The observed 
PDFs are made symmetric by taking average of the data on the left and the right hand 
sides. The measuring distances, r/rj, for the PDF of velocity fluctuations are, from the top 
to bottom in Fig. 4 of Q: 2.38, 4.76, 9.52, 19.0, 38.1, 76.2, 152, 305, 609, 1220. We see from 
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FIG. 1: Dependence of n on r/r/ 

these values that the zooming increment of the consecutive PDFs is An = — 1. The spatial 
resolution of the DNS is r m[n /r] = 2.38, and the Kolmogorov length is rj = 2.58 x 10~ 3 . 



2j, we have the 
from the top 



Adopting /i = 0.326 (q = 0.543) extracted form the analysis in Fig. 2 of 
parameters for the theoretical PDFs of velocity fluctuations in Fig. 4 of 
to bottom, (n, q') = (18.0, 1.90), (16.0, 1.80), (13.5, 1.85), (10.5, 1.75), (8.50, 1.65), ( 7.50 , 
1.60), ( 6.50 , 1.50), ( 5.50 . 1.40), ( 4.50 , 1.30), ( 3.80 , 1.20). The dependence of n on r/n is 
plotted with closed circles in Fig. E[a). The lines are adjusted by the equations [7[ 



n — —1.04 log 2 r/r) + 14.1 for r > r T , (12) 



n 



-2.391og 2 r/7? + 21.1 for r < r T . (13) 



The parameters n given above with underlines are contributing the points in Fig. da) 
adjusted by (|12jh whereas those without underlines are contributing the points in the figure 
adjusted by ([13)1 . The crossover occurs at r T /i] = 36.6, which is close to the Taylor micro 
scale A/77 = 38.3 |3] of the system. The fact that (|T2|) provides us with the correct increment 
An = — 1 indicates that the scaling exponents in Fig. 2 of 71 were extracted for r > r T with 
the interpretation that the region is the inertial range for the DNS 13 . 

Let us find out an appropriate \x value that gives the correct increment An = — 1 for the 
points represented by ()13|) . In Fig. [IJb), plotted is the dependence of n on r/n extracted 
with the appropriate fi value, fi = 0.850 (q = 0.882), for the region r < r T derived in 
accordance with the process 1) ~ 4) given above. Then, we found the parameters for the 
theoretical PDFs of velocity fluctuations with this revised fi value, from the top to bottom, 
to be (n, q') = ( 6.80 , 1.91), ( 610 , 1.90), ( 4.85 , 1.87), ( 3.72 . 1.78), ( 3.02 , 1.67), (2.53, 1.61), 
(2.15, 1.60), (1.80, 1.50), (1.50, 1.40), (1.20, 1.30). The lines given in Fig.^b) are adjusted 




by the equations 



FIG. 2: Scaling exponents 

n = —0.358 log 2 r/rj + 4.82 for r > r T , 
n = —0.995 log 2 r/rj + 8.16 for r < r T . 



(14) 
(15) 



The parameters n given with underlines are contributing the points in Fig. ^b) adjusted 
by (|15|) . whereas those without underlines are contributing the points in the figure adjusted 
by (fT4*|) . The crossover now occurs at r T /i] = 37.7. Note that PDFs with the revised /i 
value are almost the same as those given in Fig. 1 of |8] with the original fi value, i.e. it is 
impossible to see the difference in their appearance. This investigation indicates that there 
exists another range for r < r T representing a scale invariance different from the inertial 
range. A detailed investigation of the range is one of the attractive future problems. 

In Fig. |21 the transverse scaling exponents based on PDFs of velocity fluctuations ex- 



14J, Q at R A = 1054 (closed 



tracted from the time series data obtained in a wind tunnel 
circle) are analyzed by (0) with \i = 0.357 (solid line). 

In Figs. E] and |H displayed are PDFs of velocity fluctuations and of energy dissipation 
rates, respectively, extracted from the time series data ^| with the help of Taylor's 
frozen hypothesis on (a) log and (b) linear scale. The distances r/rj between two measuring 
points for PDFs of velocity fluctuations and the lengths r/rj of the region in which energy 
dissipation rates are averaged to produce its PDF are, from the top to bottom: 28.2, 56.4, 
113, 226, 451, 903. We see from these values that the increment of the consecutive PDFs is 
An = — 1. The average wind velocity in the wind tunnel is 16 m/sec. The estimated inertial 




FIG. 3: PDFs of Velocity Fluctuations 




FIG. 4: PDFs of Energy Dissipation Rates 

range is the region 44.7 < r/rj < 163. The spatial resolution is r min /r] = 14.1, and the 
Kolmogorov length is estimated as rj = 1.23 x 1CT 2 cm. For the theoretical PDFs of velocity 
fluctuations, /i = 0.380 (q = 0.610), (n, q') = (10.1, 1.78), (9.30, 1.70), (8.30, 1.65), (7.10, 

1.62) , (5.90, 1.55), (4.90, 1.50), whereas for the PDFs of energy dissipation rates, pi = 0.350 
(q = 0.574), (n, q', 9) = (7.60, 1.66, 1.81), (6.60, 1.20, 1.91), (5.60, 1.50, 1.77), (4.70, 1.68, 

1.63) , (4.20, 1.90, 1.44), (3.56, 2.30, 1.10). The values n both of velocity fluctuations and of 
energy dissipation rates have been extracted following the process 1) ~ 4) by adjusting the 
zooming increment to be An = — 1 in the analyses of PDFs. The dependence of n on r/rj 
for velocity fluctuations and for energy dissipation rates are, respectively, adjusted by the 
equations n = —1.00 log 2 r/rj + 15.0 and n = —1.00 log 2 r/rj + 12.4. Note that the scaling 
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FIG. 5: The Scaling exponents via PDF 

exponents with the derived \i values explain the extracted one in Fig. |21 within the error 
bars. 

In Fig. solid line represents the scaling exponents ( m derived by the new route via 
observed PDFs of velocity fluctuations with the help of (jOJ). The left hand side of (jHJ) is 
calculated with observed PDF data made up the lack of data for larger x n by the theoretical 
PDF n^g(a; n ). Note that, within A&A model, the difference between IT^ (x n ) and Il^ n g(x n ) 
for x n > x* n is neglected, where x* n is the point corresponding to the connection point £*. 
In the calculation of 7^, only observed data is used for 11^ (a; n ), since the integral in (JTUJl 
stops at the upper limit x* n which is about the order of the standard deviation. The formula 
for the scaling exponents (j2J) with fi = 0.380 is shown in Fig. El by dotted line which almost 
overlaps with solid line. Dashed line is the result obtained by the formula (jUJ) without making 
up the lack of data for larger x n by the substitution of theoretical PDF Il^ n g(x n ). 

The PDFs of energy transfer rates and of energy dissipation rates measured in the DNS 
on 4096 3 mesh size by Kaneda's group |3| at Ra = 1132 are, successfully, analyzed by the 
present theoretical PDFs (JOJ) and ((Zj) with = 3 for energy transfer rates and by PDFs 
(JOJ) and (jHJ with = 3 for energy dissipation rates. The observed PDFs of energy transfer 
rates are made symmetric by averaging the data on the left and the right hand sides. The 
measuring distances, r/rj, for the PDFs both of transfer rates and of dissipation rates are 
13.7, 78.1, 449. The inertial range is estimated as 63 < r/r] < 224, the Kolmogorov length 
as i] = 5.12 x 10 -4 , and the Taylor micro scale as X/i] = 66.2 jl6j. For the theoretical 
PDFs of energy transfer rates, n = 0.320 (q = 0.534), (ra, q') = (9.00, 1.75), (6.50, 1.70), 
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(3.80, 1.50), whereas for the PDFs of energy dissipation rates, \i = 0.350 (q = 0.574), (n, 
q', 6) = (7.20, 1.60, 1.30), (4.80, 1.10, 1.70), (2.30, 1.10, 4.00). The dependence of n on 
r/rj for energy transfer rates and energy dissipation rates are, respectively, adjusted by the 
equations n = — 1.04 log 2 r/rj + 12.9 and n = — 1.01 log 2 r/rj + 11.2. Note that the value of 
H has been extracted by adjusting the zooming increment to be An = —2.5 in the analyses 
of PDFs. 

IV. CONCLUSIONS 

It is shown that the formulae of PDFs within A&A model for velocity fluctuations, 
energy dissipation rates and energy transfer rates explain, successfully, corresponding PDFs 
extracted from DNSs on 2048 3 and 4096 3 mesh sizes, and experiments performed in a wind 
tunnel. It is also shown that the formulae of PDFs with appropriate zooming increment An 
corresponding to experimental situation give a new route to obtain the scaling exponents of 
velocity structure function, including intermittency exponent. 

Extracting ( 3q - for — oo < q < oo through the formula (|x n | ? ), given by @ with k = 1, out 
of the PDFs of energy dissipation rates (0 = 3), we can obtain the generalized dimension 
Dg through (|TT|) . It may be a direct proof of the multifractal distribution of singularities in 
real space. Results will be given elsewhere. Note that, for q < 0, contribution of the first 
term in becomes conspicuous, i.e., we should subtract the contributions originated from 
the term violating the invariance of the scale transformation. 

The authors would like to thank Dr. K. Yoshida for enlightening discussions. They are 
also grateful to Dr. H. Mouri and to Profs. Y. Kaneda and T. Ishihara for their kindness to 
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